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Abstract 

This work examines the effects of magnetic frustration due to competing ferromagnetic and 
antiferromagnetic Heisenberg interactions on the spin dynamics of the double-exchange model. 
When the local moments are non-colinear, a charge-density wave forms because the electrons prefer 
to sit on lines of sites that are coupled ferromagnetically. With increasing hopping energy, the local 
spins become aligned and the average spin-wave stiffness increases. Phase separation is found only 
within a narrow range of hopping energies. Results of this work are applied to the field- induced 
jump in the spin-wave stiffness observed in the manganite Pri-^^CaajMnOs with 0.3 < x < 0.4. 
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The persistence of antiferromagnetic (AFM) short-range order below the Curie temper- 
arture Tq of the manganites has been known for many years ^. Close to but below Tq, 
metallic manganites like Lao.rCao.sMnOa contain polaronic regions j^S] that are responsi- 
ble for the coexistence of propagating and diffusive spin dynamics 4]. In the remarkable 
compound Pri_a;Ca2;Mn03 with 0.3 < x < 0.4, the low-temperature ferromagnetic (FM) 
insulating phase was originally believed to be a canted AFM fCAF) 0, 0] but probably 
conta... .eg.ons w.th both FM and AFM short-.ange o.de. M When an apphed field 
B exceeds about 3 T, the resistivity drops by several orders of mag nitude M, the AFM 
regions shrink 1^1, and the spin- wave (SW) stiffness Dsw jumps by a factor of 3 p. Despite 
the recognition that short-range AFM order plays a central role in the manganites, little is 
known theoretically about how long-wavelength SW's are affected by propagating through 
both FM and AFM regions. Because electron hopping is hampered by the misalignmnent of 
the local moments 10|, AFM interactions may be expected to suppress the contribution of 
electron-mediated double-exchange (DE) to the SW dynamics . This paper examines the 
effects of AFM interactions and non-colinearity on the SW dynamics of electrons coupled to 
the local moments of a generalized Villain model Our results strongly suggest 

that the jump in the SW stiffness observed jol] in Pro.yCao.aMnOa is produced by a sharp 
increase in the hopping energy at the critical field rather than by the alignment of the AFM 
regions. 

As pictured in Fig. 1(a), the local moments Sj of the generalized Villain model are coupled 
by the FM interaction J along the y direction and by either the FM interaction J or the 
AFM interaction — r^J along the x direction. A zero-temperature CAF phase is stabilized 
when 7] exceeds Vr, which is 1/3 when B = Bz = but increases as B increases. Villain's 



original model 



12j | set rj = 1, which is the condition for full frustration. Due to the different 



environments of the a and b sites, the angle 6b at the b sites is always larger than 6a at the 
a sites, as shown in Fig. 1(b). 

Within our hybrid model, the Heisenberg interactions between the local moments are 
given by the generalized Villain model while electrons with density p = 1 — x are FM 
coupled to the local moments by Hund's coupling Jh and hop between neighboring sites with 
energy t. The DEV model (so called because it combines the DE and generalized Villain 
models) provides several advantag GS clS db basis for understanding the effects of magnetic 
frustration on the spin dynamics. First, it is one of the simplest periodic models that is 
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FIG. 1: (a) The generalized Villain model with Heisenberg couplings J or — r/J, (b) the local 
moments in the xz plane subtend angles 9a and Of, with the z axis, and (c) the electron spins also 
lie in the xz plane but subtend angles ipa < Oa and tpb < with the z axis. 

magnetically frustrated, which can be controlled through the parameter rj. In contrast to 
;he case in a DE model with AFM interactions between all neighboring local moments 
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16|, a homogeneous CAP phase is stable against phase separation [17[ except in a very 



narrow region of parameter space. Second, unlike a model with AFM exchange only, the 
DEV model supports FM order even when t = and B = 0. So it can be used to track the 
change in SW stiffness Dg^ as the electrons become mobile. Third, because it contains both 
FM and AFM Heisenberg interactions, the DEV model can be used to study insulating 
manganites like Pro.66Cao.34Mn03, where the AFM interactions arise from superexchange 
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and the FM interactions from short-range orbital and polaronic order 

For simplicity, our model is translationally symmetric with the FM and AFM Heisen- 
berg couplings arranged periodically in two dimensions. In the low-temperature phase of 
Pro.erCao.ssMnOs, the FM interactions may be confined to two-dimensional sheets in a "red 
cabbage" structure But for wavelengths longer than the thickness ~ 25 A of the FM 
sheets, the SW's will average over the FM and AFM regions. So the DEV model will pro- 
vide qualitiatively accurate predictions for the average SW stiffness D^^ = {D^^ + D^^)/2, 
which is defined in the long- wavelength limit. 
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The Hamiltonian of the DEV model is 

H=-tY. E(4Sa + c]«qj - 2 Jh E • S, - ^ J,,S, -S.-B ^(5,, + (1) 

where cj^, and are the creation and destruction operators for an electron with spin a at 
site i, Sj = (l/2)4a(7a/3Cj^ is the electronic spin, and Sj is the spin of the local moment with 
magnitude 5". The Heisenberg interactions Jij take the values J or — 77J, as described in 
Fig. 1(a). This model will be solved at zero temperature to lowest order inl/S. To guarantee 
that the contributions to the SW frequencies from DE hopping t and from the Heisenberg 
interactions Jij are of the same order in l/S", t is considered to be of the same order in 
1/5* as Jh*?, JS^ and BS (although their relative values can be quite different). Thus, the 
dimensionless parameters of the DEV model are t' — t/JS"^, rj, B' — B/JS, and J-^/JS. 
To lowest order in 1/5", the magnetic field B only couples to the local moments. While 
the theory developed below can be extended to treat all values of the Hund's couphng, for 
simplicity we shall consider the limit of large J-^S or in dimensionless terms, Jh/ JS ^ 1 
and JYiS/t > 1. 

To solve this model, a Holstein-Primakoff expansion is first performed within the rotated 
reference frame of each spin; — S a\ai, Si+ = V^Stti, and Si^ = a/ZSoI- In terms of 
electronic creation and destruction operators c^^)] and c^^^^ in the rotated reference frame of 
the local moments, the zeroth-order term (in powers of 1/y/S) in the Hamiltonian can be 
written as HQ — Efi + Hi, where to lowest order in t/ J-^S, 

Eh = ]^NJS'^ I - cos 2e„ + 77 cos 2^6 - 2 cos{ea - 9b) - B' (cos 6 a + cos ^fe) I , (2) 

= Y.\ (-JnSa - 2t COS COS Oa) + 4a^4a {'JuSa - 2t COS COS Ob) 

+ ct^^cS)2tcosA;,cos((^„ - e,)/2)]. (3) 

Here, the lattice constant is set to 1 and a = ±1 corresponds to spin up or down in the local 
reference frames. 

The electronic Hamiltonian Hi, is easily transformed into the diagonal form Hi, — 

Ek,a,r ^£4^ c?k2 by the rotations 4a = Wk^4a + «k^4a and c^i = «?^4a - Wk^4a 
where e^i = —JuSa + 4"''^ = 1 ~ '^^k''^ = (1 + (cos6'a — cos6'b) coskx/w]i)/2, 

e^^ — —t cos A;j;(cos 6a + cos 61,) ^ twk, (4) 
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W]i = y {cosOa — cos6b)'^ cos^ kx + 4cos2((6'a — 9h)/2) cos^ ky. (5) 

The =F signs refer to the r = a and h bands, respectively, and the first Brillouin zone extends 
from — vr to vr in the k^ direction but only from —tt/2 to 7r/2 in the ky direction due to the 
reduced symmetry. 

In the limit of large J-^S, the zeroth-order energy Eq = Eh + {Hh) can readily be minimized 
with respect to the angles 9a and 9b. When t = B = 0, 9b = 39a for all t]. For a fixed rj and 
B', the equilibrium angles decrease with increasing t' and 9b < 39a- The phase boundary 
between the CAF and FM phases satisfies the condition 

B'-2r] + 4 + 3Eke/4 _ 2^ {1 + r]y + {1 + E^,/8JS^y = 0, (6) 

where Eke = -((4"^) + (eif^))/2 > is the average kinetic energy of the electrons in the FM 
phase. For p = 0.66, rj = 3, and B = 0, the dependence of the equilibrium angles on t' is 
plotted in Fig.2(a). Also shown is the average spin M = S{cos9a + cos9b)/2 of the local 
moments. 

Surprisingly, the electronic occupation of the a and b sites are different with most of the 
electrons sitting on the a sites. The fraction fa of such electrons, also plotted in Fig. 2 (a), has 
a maximum of 0.58 as t' ^ and approaches 1/2 as t' — ~ 21.2. This behavior is easy 
to understand: the largest angles between neighboring spins are along the x axis between b 
sites with angles differing by 29b. When an electron hops onto a b site, it cannot easily hop 
to other b sites and so quickly moves onto a neighboring a site, where it can readily travel 
between other a sites with angular difference 29a ^ 2^;,. Hence, the non-colinearity of the 
local moments quite naturally produces a charge- density wave (CDW) with a substantial 
amplitude. A CDW with the same period as the one predicted here has in fact been observed 
in the insulating phase of Pro.7Cao.3Mn03 0, Q. 

Another surprise is that phase separation occurs within a narrow range of t' around 
10.0. Phase separation is easily seen in a plot of filling p versus chemical potential /i as a 
discontinuity Ap in It appears at fixed p as iurnps in the equilibrium angles 9,. and 

electron fraction /„. Like the Pomeranchuk instability [2l| in the two-dimensional Hubbard 
model, the phase instability in the DEV model occurs close to a Van Hove filling and is 
marked by a change in Fermi surface (FS) topology from closed to open, as shown in the 
inset to Fig. 3 where the FS is sketched for values of t' on either side of the phase-separated 
range. For t' = 10.2, the extra electrons in the neck of the a FS around k = are cancelled 
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FIG. 2: (a) The angles Oa and 9}, of the local moments, the total local magnetization M/S, the 
fraction /„ of the electrons on the a sites, and (b) the SW stiffnesses for p = 0.66, rj = 3 and B = 
versus t'. 

by the holes in the b FS around k = (tt, 7r/2). However, the phase separation is extremely 
weak and for the parameters in Fig.2, Ap ^ 0.003. 

As sketched in Fig. 1(c), the equilibrium angles ipr for the electrons are not equal to the 
angles 9r of the local moments except when JnS/t — oo. For finite JuS/t, ipr < Or the 
electrons try to align their spins as much as possible. In the limit of large Hund's coupling, 
9r — ipr oc t/JuS and the electrons always exert a small torque on the local moments. 
Hence, the relationships given above for 5^2 in terms of rf^a should contain admixtures of 
opposite-spin terms like {t/ J}iS)d^^^Q_^, where Q = (tt, 0) is the AFM Bragg vector. Since 
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FIG. 3: The SW frequencies for p = 0.66, ?? = 3, 5 = 0, and various values of t' . Different energy 
scales are used on either side of k = 0. In the inset, we plot the FS for t' = 9.8 (solid) and 10.2 
(dashed), on either side of a narrow region of phase separation. 

{six) oc sin(6'i — ipi), these new terms produce a correction to the Hund's coupling — 2 JnSj ■ Sj 
that survives in the J^S oo limit. 

After diagonalizing the band Hamiltonian iff,, the full Hamiltonian can be expanded as 
a power series in H = Hq + Hi + H2 + ■ ■ where Hi contains the torque terms 

(r) 

and is linear in the boson operators . To eliminate this first-order term and to express 
the Hamiltonian in terms of the true SW operators for the total spin Sj^tot — Sj -|- Sj, we 
perform the unitary transformation H' = e~^He^ where U is constructed to satisfy 
[U,Hq] = Hi. To lowest order in 1/5*, a trace over the Fermion degrees of freedom yields 
the modified second-order Hamiltonian H2 = H2 + [U, Hi]/2 for the SW operators only: 



(7) 



It is then straightforward to diagonalize H2 to obtain the mode frequencies Wk- The torque 
terms discussed above are required to preserve rotational symmetry and the relations co'k=o = 
B and ujq = in the CAF phase. 

In the FM phase, the SW frequency is given by the analytic result 

cjk = B + JS{3-ri+ {rj - l)cosA;^) + Ekc(2 - cosA;^)/4S' 



JSMl + ?7)2(1 - cos + 4(1 + E^jSJS'^f cos2 ky 



(8) 
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An AFM component develops when loq = 0, which yields the same condition for the phase 
boundary as Eq.®. For the FM phase, the SW stiffness obtained from the long- wavelength 
expansion cuk ~ B+D^^kl+D'^^ky is simply the sum of the DE and Heisenberg contributions: 
Aw = E^e/8S + JS{1 - r])/2 and = E],J8S + JS. For J = 0, these resuhs agree with 
the SW frequencies of the DE model first obtained by Furukawa [22] . 

In the CAE phase, the SW frequency and stiffness must be solved numerically. When 
t = 0, our results agree with Saslow and Erwin [l^ for the generalized Villain model. Results 
for are plotted in Fig. 3 for p = 0.66, t] = 3, B' = 0, and various values of t'. Above 
the phase separation region around t' ^ 10.0 but below t'^ ^ 21.2, develops kinks that 
correspond to transitions across the neck of the a FS (Q — k 0.147rx for t' = 10.2) and 
the length of the 6 FS (Q - k ^ 0.3l7rx for t' = 10.2). 

As plotted in Fig.2(b), the SW stiffness in the x direction reaches a minimum at t'^, above 
which both D^^ and D^^ are linearly increasing functions of t'. The stiffnesses in the x and 
y directions cross in the region of phase separation, where the SW's are isotropic in the 
long- wavelength limit. Notice that D^^ increases by roughly a factor of 2 as t' increases 
from zero to t'^ and the system transforms from a CAF with a CDW into a FM. 

By contrast, the effect of a magnetic field is quite different. After a sudden increase of 
the SW stiffness for very small fields that occurs in any CAE 13, 

there is a gradual increase 

in D^^ and decrease in D^^ as B increases to B^. For t' = 3 and 1] = 2, drops from 
1.275* at i? = to 0.5AJS at Be- A field also very quickly eliminates the region of phase 
separation. 



These results clearly indicate that the jump in SW stiffness observed [3] in PrcrCacsMnOa 
at a field of 3 T cannot be produced by simply aligning the AFM regions while keeping the 
bandwidth ~ t fixed. The identical resistivities in the metal-insulator transition produced by 



either a magnetic field or x-rays 



2J| suggest a common mechanism: the excitation of charge 



av 

SW 



carriers out of polaronic traps formed by the electron-lattice coupling. The doubling of D, 
found in Fig.2(b) provides strong support for this scenario. Since the integrated optical 
weight is proportional to -Eke ~ t, the jump in the hopping energy at 3 T should be reflected 
in the optical conductivity. Measurements by Okimoto et al. on Pro.6Cao.4Mn03 do reveal 
a large increase in (j{uj) and a rapid dro p i n the CDW gaps near the critical field. If the 



percolation threshold for the FM regions 



25l | is exceeded when the electrons delocalize, then 



the jump in the SW stiffness will coincide with the metal-insulator transition. Otherwise 
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the metal-insulator transition will occur at a slightly higher field. 

Because it requires two sublattices with filling x = 0.5, local CE-type AFM ordering 
in the manganites Pro.6Cao.4Mn03 j6|,|20| and Lao.rCao.sMnOa would be simplified if 

the polaronic regions were rich in holes and poor in electrons. The DEV model provides a 
natural explanation for this behavior, since the electronic fraction on b sites is substantially 
smaller than the fraction on a sites as the electrons avoid regions with more pronounced 
AFM order. 

To conclude, we have studied the effect of AFM interactions on the magnetic order and 
SW dynamics of electrons interacting with the local moments of a generalized Villain model. 
This model contains rich physics that provides insight into the SW dynamics of any itinerant 
system with competing FM and AFM Heisenberg interactions. 
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